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ABSTRACT 

High-resolution, CO(l— >0) interferometry of the ultraluminous infrared galaxy IRAS 14348-1447 is 
presented. The merger system has a molecular gas mass of ~ 3 x 10 10 M Q and a projected nuclear 
separation of 4.8 kpc (3.5"), making it one of the most molecular gas-rich galaxies known and an ideal 
candidate for studying the intermediate stages of an ultraluminous merger event. The CO morphology 
shows two molecular gas components associated with the stellar nuclei of the progenitors, consistent 
with the idea that the molecular disks are gravitationally bound by the dense bulges of the progenitor 
galaxies as the interaction proceeds. In contrast, less luminous infrared galaxies observed to date with 
projected nuclear separations of <; 5 kpc show a dominant CO component between the stellar nuclei. 
This discrepancy may be an indication that the progenitors of mergers with lower infrared luminosity do 
not possess massive bulges, and that the gas is stripped during the initial encounter of their progenitors. 
A comparison of the CO and radio luminosities of the NE and SW component show them to have 
comparable radio and CO flux ratios of f(NE)/f(SW) ~ 0.6, possibly indicating that the amount of 
star-forming molecular gas in the progenitors is correlated with the supernovae rate. The estimate of 
molecular gas masses of the nuclei and the extent of the radio emission are used to infer that the nuclei 
of IR 14348-1447 have gas densities comparable to the cores of elliptical galaxies. 

Subject headings: galaxies: ISM — infrared: galaxies — ISM: molecules — radio lines: galaxies — galaxies: 
active— individual: IRAS 14348-1447 



1. INTRODUCTION 

Galaxy tidal interactions and mergers are responsible for 
the most luminous galaxy phenomenon in the universe, be 
they star bursts or active galactic nuclei (AGN). The en- 
ergy sources are fueled by molecular gas, which is subject 
to gravitational torques, dynamical friction and dissipa- 
tion during the interaction. As a result, the distribution 
of young stars in the merger and the efficiency at which 
an AGN can be built or fueled depends not only on the 
amount of fuel available, but also on the Hubble-types 
of the progenitors (e.g., Mihos & Hernquist 1996; Mihos 
1999). 

The most energetic examples of mergers known are 
the ultraluminous infrared galaxies (ULIGs: £ir[8 — 
1000/im] > 10 12 L©). ULIGs have been of particular 
interest in the last 15 years, both because of their ex- 
treme starburst environment (e.g. Joseph & Wright 1985) 
and because of their possible evolutionary connection with 
QSOs and radio galaxies (Sanders et al. 1988a,b; Mirabel, 
Sanders, & Kazes 1989). 

In order to determine the properties of the progenitors 
of ULIGs and the mechanisms at work as the progeni- 
tors come under the gravitational influence of each other, 
a program has been initiated to map the distribution of 
CO(l — > 0) in several ULIGs for which the nuclei of the 
progenitors have yet to coalesce. In this Letter, high- 
resolution CO(l — > 0) observations of the ULIG IRAS 
14348-1447, obtained with the Owens Valley Millimeter 
Array (OVRO), are presented. With an infrared lumi- 
nosity Lir = 1.8 x 10 12 L Q and a molecular gas mass of 



~ 4.2 x 10 10 Mq (Sanders, Scoville, & Soifer 1991), IRAS 
14348-1447 is one of the most luminous and molecular gas 
rich ULIGs that show no definitive evidence of AGN ac- 
tivity (i.e, broad line emission or strong high ionization 
lines; see Veilleux, Sanders, & Kim 1997). The nuclei 
have a projected separation of ~ 3.5" (4.8 kpc: Carico 
et al. 1990; Surace, Sanders & Evans 1999; Scoville et 
al. 1999), making it ideally suited for the resolution of 
OVRO. The data presented here show evidence contrary 
to observations which have led to suggestions that molec- 
ular gas in luminous mergers as a class is stripped and 
collects between the merging stellar nuclei, but supports 
recent models that show the molecular disks in massive 
bulge galaxy-galaxy mergers are gravitationally stabilized 
against stripping (Mihos & Hernquist 1996; Mihos 1999). 
An Ho = 75 km s _1 Mpc -1 and qo = 0.5 is assumed 
throughout such that 1" subtends ~ 1.4 kpc at the red- 
shift of the galaxy (z = 0.0825). 

1.1. Interferometric Observations 

Aperture synthesis CO(l — > 0) maps of IRAS 14348- 
1447 were made with the Owens Valley Radio Observa- 
tory (OVRO) Millimeter Array during two observing pe- 
riods from 1999 February to 1999 April. The array con- 
sists of six 10.4 m telescopes, and the longest observed 
baseline was 242 m. Each telescope was configured with 
120 x 4 MHz digital correlators. During the observations, 
the nearby quasar [HB89] 1334-127 (5.89 Jy at 107 GHz; 
B1950.0 coordinates 13:34:59.81 -12:42:09.9) was observed 
every 25 minutes to monitor phase and gain variations, 
and 3C 273 and 3C 345 were observed to determine the 
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passband structure. Finally, flux calibration observations 
of Neptune were obtained. 

The OVRO data were reduced and calibrated using 
the standard Owens Valley data reduction package MMA 
(Scovillc ct al. 1992). The data were then exported to the 
mapping program DIFMAP (Shepherd, Pearson, & Taylor 
1995), and the NRAO software package AIPs was used to 
extract spectra. 

2. RESULTS 

Figure la shows a 0.8 fxm archival image of IRAS 
14348-1447 taken with the Hubble Space Telescope (HST) 
WFPC2 instrument. The merger consists of two nearly 
face-on colliding spirals with a projected nuclear separa- 
tion of 3.4" (4.8 kpc). A prominent tidal tail is visible ex- 
tending northward from the northeastern galaxy (hereafter 
IRAS 14348-1447NE), and a much fainter countertail ex- 
tends southward from the southwestern galaxy (hereafter 
IRAS 14348-1447SW). Knots, or super star clusters, are 
visible around both nuclei, as well as along the tidal tail 
of IRAS 14348-1447NE. 

Figure lb shows the CO(1-»0) emission of IRAS 14348- 
1447 in contours superposed on a three-color HST NIC- 
MOS (Near-Infrared Camera and MultiObject Spectrom- 
eter) image of the ULIG. The CO emission consists of two 
unresolved components (FWHM = 2.8" x 1.9" = 3.9 kpc 
x 2.7 kpc) which are centered on the respective stellar nu- 
clei of the progenitors. The upper limit of the CO extent 
in each galaxy is significantly less than the average effec- 
tive diameter 4 of 11.5±7.5 kpc determined from a sample 
of ~ 140 nearby spiral galaxies (Young et al. 1995). The 
total CO flux density of IRAS 14348-1447 is 40.8 Jy km 
s _1 and the CO luminosity is L' co = 8.0 x 10 9 K km s _1 
pc 2 , with ~ 42% and 68% of the emission emanating from 
the NE and SW components, respectively. The total CO 
flux density is 20% less than the flux density derived from 
the single-dish measurement of IRAS 14348-1447 with the 
NRAO 12m Telescope (half power beam width, HPBW = 
55": Sanders, Scoville, & Soifer 1991). Such a discrep- 
ancy is likely due to flux calibration uncertainties associ- 
ated with the individual observations, which can each be 
as large as 15%. Assuming a standard ratio (a) of CO lu- 
minosity to H 2 mass of 4 M Q (K km s -1 pc 2 ) -1 , which is 
similar to the value determined for the bulk of the molec- 
ular gas in the Milky Way, the total molecular gas mass 
of the pair is calculated to be 3.2 x 10 10 [a/4] M Q , or 14 
times the molecular gas mass of the Milky Way. 5 

Figure lb also shows the CO(l— >0) emission-line spec- 
tra extracted from both progenitors. The line profile of 
IRAS 14348-144 7NE appears asymmetric and blueshifted, 
consistent with the Ha line profile determined from Fabry- 
Perot observations (Mihos & Bothun 1998), and has a 
full width at half the maximum intensity line width of 
Aiifwhm ~ 350 km s" 1 . IRAS 14348-1447SW has a more 
symmetric CO line profile with a Aufwhm ~ 300 km s _1 . 
The relative velocity offset of the CO line centroids is ap- 

4 Young ct al. (1995) define the effective diameter to be the diameter that encloses 70% of the total CO emission. 

5 Radford, Solomon, &z Dowries (1991) have used theoretical models to determine that a ranges from 2—5 Mq (K km s _1 pc 2 ) -1 for a 
reasonable range of temperatures and densities. Downes & Solomon (1998) have modeled interfcromctric CO data of a sample of infrared 
luminous galaxies to derive an a = 0.8. Thus, the molecular gas mass of IRAS 14348-1447 may be a low as 6.4 X 10 9 Mq. However, note that 
the dynamical mass derived using Aiifwhm = 350 and 300 km s _1 , r = 108 and 84 pc (see §3), and assuming the reasonable inclination angle 
of the disk axis relative to the line-of-sight of 29° and 18° for the IRAS 14348-1447NE and 14348-1447SW disks, respectively, yields masses 
consistent with the molecular gas mass calculated with a = 4 Mq (K km s _1 pc 2 ) -1 . 



proximately 120 km s _1 , consistent with the offsets mea- 
sured from the optical emission lines (Veilleux et al. 1995; 
Mihos & Bothun 1998). Table 1 summarizes the properties 
discussed above. 

3. DISCUSSION 

Figure 1 clearly shows that IRAS 14348-1447 consists of 
two gas-rich spiral galaxies that have undergone gravita- 
tional response to at least one initial close approach. It 
is also clear that the molecular gas is associated with the 
individual progenitor nuclei at this stage of the merger 
process. Further, the LINER (Low Ionization Nuclear 
Emission-line Regions) emission-line spectra in both nuclei 
(Veilleux et al. 1995) indicates that shocks from either a 
low ionization AGN or supernovae resulting from massive 
starbursts with a total energy output of 1.8xl0 12 L Q have 
commenced. All of these properties can be understood in 
terms of a merger in which the molecular disks associated 
with each progenitor have been stabilized against stripping 
by a dense, massive stellar bulge. Given the upper limit 
on the gas distribution of IRAS 14348-1447 relative to lo- 
cal spiral galaxies, the interaction has likely already driven 
the molecular gas inwards, resulting in enhanced nuclear 
activity (see simulations by Mihos & Hernquist 1996 and 
Mihos 1999). 

Disk stabilization against stripping is also applicable to 
other intermediate stage ULIGs. Recent observations of 
Arp 220 have shown that the major CO concentrations 
in this merger are on the individual stellar nuclei, which 
have a projected separation of 0.95" (~ 350 pc; Sakamoto 
et al. 1999). Likewise, PKS 1345+12, a ULIG with very 
warm, Seyfert-like infrared colors relative to Arp 220 and 
IRAS 14348-1447, shows the CO emission concentrated 
only on the active radio nucleus, indicating that gas infall 
and AGN activity has been triggered by interactions with 
the companion galaxy (Evans et al. 1999). The compan- 
ion galaxy has colors consistent with an elliptical galaxy 
(Surace et al. 1998), which explains its lack of detectable 
molecular gas. 

In contrast, the morphologies of all of the lower luminos- 
ity luminous infrared galaxies (LIGs: Lir = ifj 110-11 " 
Lq) with small (<J 5 kpc) projected nuclear separations 
observed to date are markedly different; Arp 244, VV 114, 
NGC 6240, and NGC 6090 have molecular gas predomi- 
nantly in one component between the two nuclei (Stanford 
et al. 1991; Yun et al. 1994; Tacconi et al. 1999; Bryant 
& Scoville 1999). Tacconi et al. (1999) have speculated 
that the molecular gas in NGC 6240 has been ram-pressure 
stripped, and that the nuclei may later sweep up gas as the 
galaxy evolves into ULIGs such as Arp 220. An alternative 
explanation may be that the processes affecting the gas in 
LIGs differ from ULIGs because the progenitors of each 
luminosity class differ. Specifically, ULIGs may primarily 
be examples of equal mass galaxy mergers with dense stel- 
lar bulges, thus the gas is driven into the nuclear region 
of the respective progenitors as the merger advances. In 
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contrast, LIGs may be predominantly collisions of galaxies 
of different masses and relatively low-density bulges (e.g., 
note the differences in the near-infrared stellar morpholo- 
gies of the progenitor galaxies of NGC 6090 and VV 114: 
Dinshaw et al. 1999; Scoville et al. 1999), thus the gas 
is stripped from the progenitors with extreme efficiency. 
As a result, many LIGs may never evolve into ULIGs. A 
larger survey of double nuclei ULIGs, consisting of those 
with obvious AGN signatures and those without, is un- 
derway to investigate the ubiquity of these results. Such a 
survey will also benefit from kinematic determinations of 
the stellar bulge masses of the progenitors. 

A comparison of the molecular gas and radio fluxes and 
morphologies of the progenitors of IRAS 14348-1447 can 
be used to derive the gas densities of the nuclei. The mea- 
sured IRAS 14348-1447NE to IRAS 14348-1447SW flux 
density ratios of CO and the 8.44 GHz and 1.49 GHz ra- 
dio emission (i.e., Condon et al. 1990; 1991) yield values of 
f(NE)/f(SW) ~ 0.61, 0.62, and 0.58, respectively. The 
implication of this result is that the molecular gas mass of 
each component is related to the source of the radio emis- 
sion. This can be understood if the nuclear radio emis- 
sion is due to synchrotron emission from supernovae and 
if both galaxies have similar initial mass functions. Thus, 
the CO and radio flux density ratios of the progenitors are 
similar because the same fraction of massive stars are pro- 
duced per unit of star-forming molecular gas. Therefore, if 
the likely assumption is made that the extent of the radio 
emission of each nucleus (FWHM(NE) - 0.16" [220 pc] 
and FWHM(SW) - 0.12 [170 pc]) is similar to the true 
extent of the molecular gas (i.e., the supernovae fill the 
same volume as the gas they are formed from), then the 
northeastern and southwestern nuclei have gas densities of 
- 2.4 x 10 3 [a/A] and 7.7 x 10 3 [a/A] M pc~ 3 , respec- 
tively. Given the uncertainty in the value of a, the density 
and the velocity dispersion of the gas in each nucleus are 
comparable to the stellar densities and velocity dispersions 
of elliptical galaxy cores (Faber et al. 1997), supporting 
the likely connection between ULIGs and the formation of 
elliptical galaxy bulges (Kormendy & Sanders 1992). 



Molecular gas-rich ULIGs in the local universe such as 
IRAS 14348-1447 may provide insights to the nature of 
massive galaxy formation in the universe. Figure 2 shows 
a plot of the logarithm of CO (1^0) luminosity of LIGs and 
ULIGs versus their rcdshift. On this plot, IR 14348-1447 
is shown as an asterix enclosed in a circle. The rise in L' co 
at z < 0.04 is simply due to the space density of the flux- 
limited infrared luminous galaxy sample. However, the 
leveling off of L' co beyond z > 0.07 is a possible indica- 
tion that, due to self regulating processes, galaxies do not 
contain molecular gas masses in excess of ~ 4 x 10 10 [a/4] 
M Q (Evans et al. 1996; Frayer et al. 1999). The observed 
flatness of L' co beyond z ~ 0.07 remains constant out to 
rcdshifts of 4.7 (Frayer et al. 1999). Thus, in terms of the 
richness of the interstellar medium, galaxies such as IRAS 
14348-1447 appear to be the low redshift counterparts of 
molecular gas-rich, high-redshift galaxies detected over the 
last decade (e.g. see Frayer et al. 1999 for a summary). 
If a substantial fraction of these systems have nuclear gas 
densities comparable to IRAS 14348-1447, then they as a 
class are the likely progenitors of massive elliptical galax- 
ies. 
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Figure Captions 



Figure l.a) Hubble Space Telescope WFPC2 image of IRAS 14348-1447. b) CO(l — > 0) contours of the merger superim- 
posed on a three-color composite NICMOS image (bluc=l.l /im, green =1.6 /im, red=2.2 /zm). The CO contours are 
plotted as 20%, 30%, 40%, 50%, 60%, 70%, 80%, 90%, and 99% the peak flux of 0.0283 Jy/beam. The CO emission for 
each progenitor is unresolved, with a beam FWHM of 2.8"xl.9" at a position angle of -21.6°. Extracted spectra of the 
SW and NE progenitors are also shown. For the images, north is up and east is to the left. 

Figure 2. The plot of log(X' co ) versus redshift (z) for a flux-limited sample (/60/xm > 5.24 Jy) of infrared luminous galaxies 
and a sample of ultraluminous infrared galaxies. The data have been obtained from Sanders, Scoville, & Soifer (1991) 
and Solomon et al. (1997). The data point representing IRAS 14348-1447 is encircled. 
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Table 1 

Molecular Gas Properties of IRAS 14348-1447 



Component 


RA Dec 


z 


AwpWHM 


ScoAv 


M(R 2 ) b 




(B1950.0) a 




km s _1 


Jy km s -1 


xlO 10 M 


14348-1447NE 


14:34:53.42 -14:47:234 


0.0823 


350 


15.6 


1.2 


14348-1447SW 


14:34:53.30 -14:47:26.1 


0.0827 


300 


25.2 


2.0 


Total 








40.8 


3.2 



a The coordinates are taken from Condon et al. (1991). 
b The H2 masses are derived assuming a = 4 (see text). 
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